Novel signals and polygenic score for height are associated with pubertal growth traits in Southwestern American Indians.

Most genetic variants associated with adult height have been identified through large genome-wide association studies (GWASs) in European-ancestry cohorts. However, it is unclear how these variants influence linear growth during adolescence. This study uses anthropometric and genotypic data from a longitudinal study conducted in an American Indian community in Arizona between 1965-2007. Growth parameters (i.e. height, velocity, and timing of growth spurt) were derived from the Preece-Baines growth model, a parametric growth curve fitted to longitudinal height data, in 787 participants with height measurements spanning the whole period of growth. Heritability estimates suggested that genetic factors could explain 25% to 71% of the variance of pubertal growth traits. We performed a GWAS of growth parameters, testing their associations with 5 077 595 imputed or directly genotyped variants. Six variants associated with height at peak velocity (P < 5 × 10-8, adjusted for sex, birth year and principal components). Implicated genes include NUDT3, previously associated with adult height, and PACSIN1. Two novel variants associated with duration of growth spurt (P < 5 × 10-8) in LOC105375344, an uncharacterized gene with unknown function. We finally examined the association of growth parameters with a polygenic score for height derived from 9557 single nucleotide polymorphisms (SNPs) identified in the GIANT meta-analysis for which genotypic data were available for the American Indian study population. Height polygenic score was correlated with the magnitude and velocity of height growth that occurred before and at the peak of the adolescent growth spurt, indicating overlapping genetic architecture, with no influence on the timing of adolescent growth.

Epistasis Between HLA-DRB1*16:02:01 and SLC16A11 T-C-G-T-T Reduces Odds for Type 2 Diabetes in Southwest American Indians.

We sought to identify genetic/immunologic contributors of type 2 diabetes in an indigenous American community by genotyping all study participants for both high resolution HLA-DRB1 alleles and SLC16A11 to test their risk and/or protection for T2D. These genes were selected based on independent reports that HLA-DRB1*16:02:01 is protective for T2D and that SLC16A11 associates with T2D in individuals with BMI<35kg/m2, and here test the interaction of the two loci with a more complete dataset and perform a BMI sensitivity test. We define the risk-protection haplotype of SLC16A11, T-C-G-T-T, as allele "2" of a di-allelic genetic model with three genotypes, SLC16A11*11, *12, and *22, where allele "1" is the wildtype. Both earlier findings were confirmed. Together in the same logistic model with BMI≥35, DRB1*16:02:01 remains protective, 0.73, while SLC16A11 switches from risk to protection OR = 0.57 (*22) and 0.78 (*12), respectively; an added interaction term was statistically significant (OR = 0.49 with *12). Bootstrapped (b=10,000) statistical power of interaction, 0.4801, yielded mean OR = 0.43. Sensitivity analysis demonstrated the interaction significant in BMI range 30-41. To investigate the epistasis we used the primary function of the HLA-DRB1 molecule, peptide binding and presentation, to search the entire array of 15mer peptides for both the wildtype and ancient human SLC16A11 molecules for a pattern of strong binding that was associated with risk and protection for T2D. Applying computer binding algorithms suggests the core peptide at SLC16A11 D127G, FSAFASGLL, might be key for moderating risk for T2D with potential implications for T1D.

Insulin resistance before type 2 diabetes onset is associated with increased risk of albuminuria after diabetes onset: A prospective cohort study.

AIM: Reduced renal insulin signalling is implicated in the pathogenesis of albuminuria. We sought to investigate whether insulin action and secretion, measured before diabetes onset, are associated with the development of albuminuria after diabetes onset. MATERIALS AND METHODS: Baseline body composition, insulin sensitivity by hyperinsulinaemic-euglycaemic clamp at submaximal and maximal insulin stimulation (240 and 2400 pmol/m2/min; M-low and M-high), and insulin secretion by intravenous glucose tolerance test [acute insulin response (AIR)] were measured in 170 Southwestern Indigenous American adults who subsequently developed diabetes. After diabetes onset and during the median follow-up of 13.6 years, 81 participants (48%) developed albuminuria (urine albumin-to-creatinine ratio ≥30 mg/g). Separate associations of M-low, M-high and AIR (per 1-SD change) with the risk of albuminuria were assessed by Cox regression models adjusted for age, sex and body fat (%). RESULTS: Participants who developed albuminuria were of similar age (26.4 ± 5.4 vs. 27.5 ± 6.1 years), sex (46% vs. 48% male), body fat (36.4 ± 7.5 vs. 35.7 ± 7.9%) and AIR [2.3 ± 0.3 vs. 2.3 ± 0.3, pmol/L (log)] as those who did not develop albuminuria but had lower insulin sensitivity [M-low: 0.33 ± 0.08 vs. 0.36 ± 0.12, p = .03; M-high: 0.87 ± 0.11 vs. 0.91 ± 0.12, p = .02; mg/kg-metabolic body size/min (log)]. In separate adjusted models, lower M-low and M-high were both associated with an increased risk for albuminuria [hazard ratio (HR) 1.51, 95% confidence interval (CI) 1.14, 2.00, p = .004; HR 1.31, 95% CI 1.06, 1.63, p = .01), whereas AIR was not (HR 1.15, 95% CI 0.87, 1.56, p = .3). CONCLUSIONS: Lower insulin sensitivity is associated with the development of albuminuria, suggesting a role for insulin signalling in the pathogenesis of proteinuria.

Trans-ancestral genome-wide association study of longitudinal pubertal height growth and shared heritability with adult health outcomes.

BACKGROUND: Pubertal growth patterns correlate with future health outcomes. However, the genetic mechanisms mediating growth trajectories remain largely unknown. Here, we modeled longitudinal height growth with Super-Imposition by Translation And Rotation (SITAR) growth curve analysis on ~ 56,000 trans-ancestry samples with repeated height measurements from age 5 years to adulthood. We performed genetic analysis on six phenotypes representing the magnitude, timing, and intensity of the pubertal growth spurt. To investigate the lifelong impact of genetic variants associated with pubertal growth trajectories, we performed genetic correlation analyses and phenome-wide association studies in the Penn Medicine BioBank and the UK Biobank. RESULTS: Large-scale growth modeling enables an unprecedented view of adolescent growth across contemporary and 20th-century pediatric cohorts. We identify 26 genome-wide significant loci and leverage trans-ancestry data to perform fine-mapping. Our data reveals genetic relationships between pediatric height growth and health across the life course, with different growth trajectories correlated with different outcomes. For instance, a faster tempo of pubertal growth correlates with higher bone mineral density, HOMA-IR, fasting insulin, type 2 diabetes, and lung cancer, whereas being taller at early puberty, taller across puberty, and having quicker pubertal growth were associated with higher risk for atrial fibrillation. CONCLUSION: We report novel genetic associations with the tempo of pubertal growth and find that genetic determinants of growth are correlated with reproductive, glycemic, respiratory, and cardiac traits in adulthood. These results aid in identifying specific growth trajectories impacting lifelong health and show that there may not be a single "optimal" pubertal growth pattern.

Diagnostic criteria and etiopathogenesis of type 2 diabetes and its complications: Lessons from the Pima Indians.

The Phoenix Epidemiology and Clinical Research Branch of the National Institute of Diabetes and Digestive and Kidney Diseases has conducted prospective studies of diabetes and its complications in the Pima Indians living in Arizona, USA for over 50 years. In this review we highlight areas in which these studies provided vital insights into the criteria used to diagnose type 2 diabetes, the pathophysiologic changes that accompany the development of type 2 diabetes, and the course and determinants of diabetes complications-focusing specifically on diabetic kidney disease. We include data from our longitudinal population-based study of diabetes and its complications, studies on the role of insulin resistance and insulin secretion in the pathophysiology of type 2 diabetes, and in-depth studies of diabetic kidney disease that include measures of glomerular function and research kidney biopsies. We also focus on the emerging health threat posed by youth-onset type 2 diabetes, which was first seen in the Pima Indians in the 1960s and is becoming an increasing issue worldwide.

Mechanisms of weight loss-induced remission in people with prediabetes: a post-hoc analysis of the randomised, controlled, multicentre Prediabetes Lifestyle Intervention Study (PLIS).

BACKGROUND: Remission of type 2 diabetes can occur as a result of weight loss and is characterised by liver fat and pancreas fat reduction and recovered insulin secretion. In this analysis, we aimed to investigate the mechanisms of weight loss- induced remission in people with prediabetes. METHODS: In this prespecified post-hoc analysis, weight loss-induced resolution of prediabetes in the randomised, controlled, multicentre Prediabetes Lifestyle Intervention Study (PLIS) was assessed, and the results were validated against participants from the Diabetes Prevention Program (DPP) study. For PLIS, between March 1, 2012, and Aug 31, 2016, participants were recruited from eight clinical study centres (including seven university hospitals) in Germany and randomly assigned to receive either a control intervention, a standard lifestyle intervention (ie, DPP-based intervention), or an intensified lifestyle intervention for 12 months. For DPP, participants were recruited from 23 clinical study centres in the USA between July 31, 1996, and May 18, 1999, and randomly assigned to receive either a standard lifestyle intervention, metformin, or placebo. In both PLIS and DPP, only participants who were randomly assigned to receive lifestyle intervention or placebo and who lost at least 5% of their bodyweight were included in this analysis. Responders were defined as people who returned to normal fasting plasma glucose (FPG; <5·6 mmol/L), normal glucose tolerance (<7·8 mmol/L), and HbA1c less than 39 mmol/mol after 12 months of lifestyle intervention or placebo or control intervention. Non-responders were defined as people who had FPG, 2 h glucose, or HbA1c more than these thresholds. The main outcomes for this analysis were insulin sensitivity, insulin secretion, visceral adipose tissue (VAT), and intrahepatic lipid content (IHL) and were evaluated via linear mixed models. FINDINGS: Of 1160 participants recruited to PLIS, 298 (25·7%) had weight loss of 5% or more of their bodyweight at baseline. 128 (43%) of 298 participants were responders and 170 (57%) were non-responders. Responders were younger than non-responders (mean age 55·6 years [SD 9·9] vs 60·4 years [8·6]; p<0·0001). The DPP validation cohort included 683 participants who lost at least 5% of their bodyweight at baseline. Of these, 132 (19%) were responders and 551 (81%) were non-responders. In PLIS, BMI reduction was similar between responders and non-responders (responders mean at baseline 32·4 kg/m2 [SD 5·6] to mean at 12 months 29·0 kg/m2 [4·9] vs non-responders 32·1 kg/m2 [5·9] to 29·2 kg/m2 [5·4]; p=0·86). However, whole-body insulin sensitivity increased more in responders than in non-responders (mean at baseline 291 mL/[min × m2], SD 60 to mean at 12 months 378 mL/[min × m2], 56 vs 278 mL/[min × m2], 62, to 323 mL/[min × m2], 66; p<0·0001), whereas insulin secretion did not differ within groups over time or between groups (responders mean at baseline 175 pmol/mmol [SD 64] to mean at 12 months 163·7 pmol/mmol [60·6] vs non-responders 158·0 pmol/mmol [55·6] to 154·1 pmol/mmol [56·2]; p=0·46). IHL decreased in both groups, without a difference between groups (responders mean at baseline 10·1% [SD 8·7] to mean at 12 months 3·5% [3·9] vs non-responders 10·3% [8·1] to 4·2% [4·2]; p=0·34); however, VAT decreased more in responders than in non-responders (mean at baseline 6·2 L [SD 2·9] to mean at 12 months 4·1 L [2·3] vs 5·7 L [2·3] to 4·5 L [2·2]; p=0·0003). Responders had a 73% lower risk of developing type 2 diabetes than non-responders in the 2 years after the intervention ended. INTERPRETATION: By contrast to remission of type 2 diabetes, resolution of prediabetes was characterised by an improvement in insulin sensitivity and reduced VAT. Because return to normal glucose regulation (NGR) prevents development of type 2 diabetes, we propose the concept of remission of prediabetes in analogy to type 2 diabetes. We suggest that remission of prediabetes should be the primary therapeutic aim in individuals with prediabetes. FUNDING: German Federal Ministry for Education and Research via the German Center for Diabetes Research; the Ministry of Science, Research and the Arts Baden-Württemberg; the Helmholtz Association and Helmholtz Munich; the Cluster of Excellence Controlling Microbes to Fight Infections; and the German Research Foundation.

Glycemic Measures in Childhood as Predictors of Future Diabetes-Related Microvascular Complications in an Indigenous American Population.

OBJECTIVE: To examine the role of glycemic measures performed during childhood in predicting future diabetes-related nephropathy and retinopathy in a high-risk indigenous American cohort. RESEARCH DESIGN AND METHODS: We studied associations between glycated hemoglobin (HbA1c) and 2-h plasma glucose (PG), measured during childhood (age 5 to <20 years) in a longitudinal observational study of diabetes and its complications (1965-2007), and future albuminuria (albumin creatinine ratio [ACR] ≥30 mg/g), severe albuminuria (ACR ≥300 mg/g), and retinopathy (at least one microaneurysm or hemorrhage or proliferative retinopathy on direct ophthalmoscopy). Areas under the receiver operating characteristic curve (AUCs) for childhood glycemic measures when predicting nephropathy and retinopathy were compared. RESULTS: Higher baseline levels of HbA1c and 2-h PG significantly increased the risk of future severe albuminuria (HbA1c: hazard ratio [HR] 1.45 per %; 95% CI 1.02-2.05 and 2-h PG: HR 1.21 per mmol/L; 95% CI 1.16-1.27). When categorized by baseline HbA1c, children with prediabetes had a higher incidence of albuminuria (29.7 cases per 1,000 person-years [PY]), severe albuminuria (3.8 cases per 1,000 PY), and retinopathy (7.1 cases per 1,000 PY) than children with normal HbA1c levels (23.8, 2.4, and 1.7 cases per 1,000 PY, respectively); children with diabetes at baseline had the highest incidence of the three complications. No significant differences were observed between AUCs for models with HbA1c, 2-h PG, and fasting PG when predicting albuminuria, severe albuminuria, or retinopathy. CONCLUSIONS: In this study, higher glycemia levels ascertained by HbA1c and 2-h PG during childhood were associated with future microvascular complications; this demonstrates the potential utility of screening tests performed in high-risk children in predicting long-term health outcomes.

DNA methylation markers for kidney function and progression of diabetic kidney disease.

Epigenetic markers are potential biomarkers for diabetes and related complications. Using a prospective cohort from the Hong Kong Diabetes Register, we perform two independent epigenome-wide association studies to identify methylation markers associated with baseline estimated glomerular filtration rate (eGFR) and subsequent decline in kidney function (eGFR slope), respectively, in 1,271 type 2 diabetes subjects. Here we show 40 (30 previously unidentified) and eight (all previously unidentified) CpG sites individually reach epigenome-wide significance for baseline eGFR and eGFR slope, respectively. We also develop a multisite analysis method, which selects 64 and 37 CpG sites for baseline eGFR and eGFR slope, respectively. These models are validated in an independent cohort of Native Americans with type 2 diabetes. Our identified CpG sites are near genes enriched for functional roles in kidney diseases, and some show association with renal damage. This study highlights the potential of methylation markers in risk stratification of kidney disease among type 2 diabetes individuals.

Adolescent Growth Spurt and Type 2 Diabetes Risk in Southwestern American Indians.

Early puberty onset is associated with higher risk of diabetes, but most studies have not accounted for childhood factors that may confound the association. Using data from a study conducted in an Indigenous community in Arizona (1965-2007), we examined associations of timing and velocity of the adolescent growth spurt with type 2 diabetes, and whether these associations are mediated by childhood body mass index and insulinemia. Adolescent growth parameters were derived from the Preece-Baines growth model, a parametric growth curve fitted to longitudinal height data, for 861 participants with height measurements spanning the whole period of growth. In males, older age at take-off, age at peak velocity, and age at maturation were associated with decreased prevalence of diabetes (odds ratio (OR) = 0.43 per year, 95% confidence interval (CI): 0.27, 0.69; OR = 0.50, 95% CI: 0.35, 0.72; OR = 0.58, 95% CI: 0.41, 0.83, respectively), while higher velocity at take-off was associated with increased risk (OR = 3.47 per cm/year, 95% CI: 1.87, 6.42) adjusting for age, birth year, and maternal diabetes. Similar results were observed with incident diabetes. Our findings suggest that an early and accelerated adolescent growth spurt is a risk factor for diabetes, at least in males. These associations are only partially explained by measures of adiposity and insulinemia.

The utility of a type 2 diabetes polygenic score in addition to clinical variables for prediction of type 2 diabetes incidence in birth, youth and adult cohorts in an Indigenous study population.

AIMS/HYPOTHESIS: There is limited information on how polygenic scores (PSs), based on variants from genome-wide association studies (GWASs) of type 2 diabetes, add to clinical variables in predicting type 2 diabetes incidence, particularly in non-European-ancestry populations. METHODS: For participants in a longitudinal study in an Indigenous population from the Southwestern USA with high type 2 diabetes prevalence, we analysed ten constructions of PS using publicly available GWAS summary statistics. Type 2 diabetes incidence was examined in three cohorts of individuals without diabetes at baseline. The adult cohort, 2333 participants followed from age ≥20 years, had 640 type 2 diabetes cases. The youth cohort included 2229 participants followed from age 5-19 years (228 cases). The birth cohort included 2894 participants followed from birth (438 cases). We assessed contributions of PSs and clinical variables in predicting type 2 diabetes incidence. RESULTS: Of the ten PS constructions, a PS using 293 genome-wide significant variants from a large type 2 diabetes GWAS meta-analysis in European-ancestry populations performed best. In the adult cohort, the AUC of the receiver operating characteristic curve for clinical variables for prediction of incident type 2 diabetes was 0.728; with the PS, 0.735. The PS's HR was 1.27 per SD (p=1.6 × 10-8; 95% CI 1.17, 1.38). In youth, corresponding AUCs were 0.805 and 0.812, with HR 1.49 (p=4.3 × 10-8; 95% CI 1.29, 1.72). In the birth cohort, AUCs were 0.614 and 0.685, with HR 1.48 (p=2.8 × 10-16; 95% CI 1.35, 1.63). To further assess the potential impact of including PS for assessing individual risk, net reclassification improvement (NRI) was calculated: NRI for the PS was 0.270, 0.268 and 0.362 for adult, youth and birth cohorts, respectively. For comparison, NRI for HbA1c was 0.267 and 0.173 for adult and youth cohorts, respectively. In decision curve analyses across all cohorts, the net benefit of including the PS in addition to clinical variables was most pronounced at moderately stringent threshold probability values for instituting a preventive intervention. CONCLUSIONS/INTERPRETATION: This study demonstrates that a European-derived PS contributes significantly to prediction of type 2 diabetes incidence in addition to information provided by clinical variables in this Indigenous study population. Discriminatory power of the PS was similar to that of other commonly measured clinical variables (e.g. HbA1c). Including type 2 diabetes PS in addition to clinical variables may be clinically beneficial for identifying individuals at higher risk for the disease, especially at younger ages.

Identification and characterization of the long non-coding RNA NFIA-AS2 as a novel locus for body mass index in American Indians.

BACKGROUND: Genome-wide association studies have shown that body mass index (BMI), an estimate of obesity, is highly polygenic. Individual variants typically have small effect sizes, making it challenging to identify unique loci in under-represented ethnic groups which lack statistical power due to their small sample size. Yet obesity is a major health disparity and is particularly prevalent in southwestern American Indians. Here, we identify and characterize a new locus for BMI that was detected by analyzing moderate associations with BMI obtained in a population-based sample of southwestern American Indians together with the well-powered GIANT dataset. METHODS: Genotypes for 10.5 million variants were tested for association with BMI in 5870 American Indians and 2600 variants that showed an association P < 10-3 in the American Indian sample were combined in a meta-analysis with the BMI data reported in GIANT (N = 240,608). The newly identified gene, NFIA-AS2 was functionally characterized, and the impact of its lead associated variant rs1777538 was studied both in-silico and in-vitro. RESULTS: Rs1777538 (T/C; C allele frequency = 0.16 in American Indians and 0.04 in GIANT, meta-analysis P = 5.0 × 10-7) exhibited a large effect in American Indians (1 kg/m2 decrease in BMI per copy of C allele). NFIA-AS2 was found to be a nuclear localized long non-coding RNA expressed in tissues pertinent to human obesity. Analysis of this variant in human brown preadipocytes showed that NFIA-AS2 transcripts carrying the C allele had increased RNA degradation compared to the T allele transcripts (half-lives = 9 h, 13 h respectively). During brown adipogenesis, NFIA-AS2 featured a stage-specific regulation of nearby gene expression where rs1777538 demonstrated an allelic difference in regulation in the mature adipocytes (the strongest difference was observed for L1TD1, P = 0.007). CONCLUSION: Our findings support a role for NFIA-AS2 in regulating pathways that impact BMI.

Identification of Genetic Variation Influencing Metformin Response in a Multiancestry Genome-Wide Association Study in the Diabetes Prevention Program (DPP).

Genome-wide significant loci for metformin response in type 2 diabetes reported elsewhere have not been replicated in the Diabetes Prevention Program (DPP). To assess pharmacogenetic interactions in prediabetes, we conducted a genome-wide association study (GWAS) in the DPP. Cox proportional hazards models tested associations with diabetes incidence in the metformin (MET; n = 876) and placebo (PBO; n = 887) arms. Multiple linear regression assessed association with 1-year change in metformin-related quantitative traits, adjusted for baseline trait, age, sex, and 10 ancestry principal components. We tested for gene-by-treatment interaction. No significant associations emerged for diabetes incidence. We identified four genome-wide significant variants after correcting for correlated traits (P < 9 × 10-9). In the MET arm, rs144322333 near ENOSF1 (minor allele frequency [MAF]AFR = 0.07; MAFEUR = 0.002) was associated with an increase in percentage of glycated hemoglobin (per minor allele, ß = 0.39 [95% CI 0.28, 0.50]; P = 2.8 × 10-12). rs145591055 near OMSR (MAF = 0.10 in American Indians) was associated with weight loss (kilograms) (per G allele, ß = -7.55 [95% CI -9.88, -5.22]; P = 3.2 × 10-10) in the MET arm. Neither variant was significant in PBO; gene-by-treatment interaction was significant for both variants [P(G×T) < 1.0 × 10-4]. Replication in individuals with diabetes did not yield significant findings. A GWAS for metformin response in prediabetes revealed novel ethnic-specific associations that require further investigation but may have implications for tailored therapy.

Effects of the ABCC8 R1420H loss-of-function variant on beta-cell function, diabetes incidence, and retinopathy.

INTRODUCTION: The ABCC8 gene regulates insulin secretion and plays a critical role in glucose homeostasis. The effects of an ABCC8 R1420H loss-of-function variant on beta-cell function, incidence of type 2 diabetes, and age-at-onset, prevalence, and progression of diabetes complications were assessed in a longitudinal study in American Indians. RESEARCH DESIGN AND METHODS: We analyzed beta-cell function through the relationship between insulin secretion and insulin sensitivity in members of this population without diabetes aged ≥5 years using standard major axis regression. We used hierarchical logistic regression models to study cross-sectional associations with diabetes complications including increased albuminuria (albumin-to-creatinine ratio (ACR) ≥30 mg/g), severe albuminuria (ACR ≥300 mg/g), reduced estimated glomerular filtration rate (eGFR <60 mL/min/1.73 m2), and retinopathy. This study included 7675 individuals (254 variant carriers) previously genotyped for the R1420H with available phenotypic data and with a median follow-up time of 13.5 years (IQR 4.5-26.8). RESULTS: Variant carriers had worse beta-cell function than non-carriers (p=0.0004; on average estimated secretion was 22% lower, in carriers), in children and adults, with no difference in insulin sensitivity (p=0.50). At any body mass index and age before 35 years, carriers had higher type 2 diabetes incidence. This variant did not associate with prevalence of increased albuminuria (OR 0.87, 95% CI 0.66 to 1.16), severe albuminuria (OR 0.96, 95% CI 0.55 to 1.68), or reduced eGFR (OR 0.44, 95% CI 0.18 to 1.06). By contrast, the variant significantly associated with higher retinopathy prevalence (OR 1.74, 95% CI 1.19 to 2.53) and this association was only partially mediated (<11%) by glycemia, duration of diabetes, risk factors of retinopathy, or insulin use. Retinopathy prevalence in carriers was higher regardless of diabetes presence. CONCLUSIONS: The ABCC8 R1420H variant is associated with increased risks of diabetes and of retinopathy, which may be partially explained by higher glycemia levels and worse beta-cell function.

Weight Loss, Lifestyle Intervention, and Metformin Affect Longitudinal Relationship of Insulin Secretion and Sensitivity.

CONTEXT: Insulin secretion and sensitivity regulate glycemia, with inadequately compensated deficiencies leading to diabetes. OBJECTIVE: We investigated effects of weight loss, an intensive lifestyle intervention (ILS), and metformin on the relationship between insulin secretion and sensitivity using repository data from 2931 participants in the Diabetes Prevention Program clinical trial in adults at high risk of developing type 2 diabetes. METHODS: Insulin secretion and sensitivity were estimated from insulin and glucose concentrations in fasting and 30-minute postload serum samples at baseline and 1, 2, and 3 years after randomization, during the active intervention phase. The nonlinear relationship of secretion and sensitivity was evaluated by standardized major axis regression to account for variability in both variables. Insulin secretory demand and compensatory insulin secretion were characterized by distances along and away from the regression line, respectively. RESULTS: ILS and metformin decreased secretory demand while increasing compensatory insulin secretion, with greater effects of ILS. Improvements were directly related to weight loss; decreased weight significantly reduced secretory demand (b=-0.144 SD; 95% CI (-0.162, -0.125)/5 kg loss) and increased compensatory insulin secretion (b = 0.287 SD, 95% CI (0.261, 0.314)/5 kg loss). In time-dependent hazard models, increasing compensatory insulin secretion (hazard ratio [HR] = 0.166 per baseline SD, 95% CI 0.133, 0.206) and weight loss (HR = 0.710 per 5 kg loss, 95% CI 0.613, 0.819) predicted lower diabetes risk. CONCLUSION: Diabetes risk reduction was directly related to the amount of weight loss, an effect mediated by lowered insulin secretory demand (due to increased insulin sensitivity) coupled with improved compensatory insulin secretion.

Association of protein function-altering variants with cardiometabolic traits: the strong heart study.

Clinical and biomarker phenotypic associations for carriers of protein function-altering variants may help to elucidate gene function and health effects in populations. We genotyped 1127 Strong Heart Family Study participants for protein function-altering single nucleotide variants (SNV) and indels selected from a low coverage whole exome sequencing of American Indians. We tested the association of each SNV/indel with 35 cardiometabolic traits. Among 1206 variants (average minor allele count = 20, range of 1 to 1064), ~ 43% were not present in publicly available repositories. We identified seven SNV-trait significant associations including a missense SNV at ABCA10 (rs779392624, p = 8 × 10-9) associated with fasting triglycerides, which gene product is involved in macrophage lipid homeostasis. Among non-diabetic individuals, missense SNVs at four genes were associated with fasting insulin adjusted for BMI (PHIL, chr6:79,650,711, p = 2.1 × 10-6; TRPM3, rs760461668, p = 5 × 10-8; SPTY2D1, rs756851199, p = 1.6 × 10-8; and TSPO, rs566547284, p = 2.4 × 10-6). PHIL encoded protein is involved in pancreatic ß-cell proliferation and survival, and TRPM3 protein mediates calcium signaling in pancreatic ß-cells in response to glucose. A genetic risk score combining increasing insulin risk alleles of these four genes was associated with 53% (95% confidence interval 1.09, 2.15) increased odds of incident diabetes and 83% (95% confidence interval 1.35, 2.48) increased odds of impaired fasting glucose at follow-up. Our study uncovered novel gene-trait associations through the study of protein-coding variants and demonstrates the advantages of association screenings targeting diverse and high-risk populations to study variants absent in publicly available repositories.

Functional characterization of a novel p.Ser76Thr variant in IGFBP4 that associates with body mass index in American Indians.

Insulin-like growth factor binding protein 4 (IGFBP4) is involved in adipogenesis, and IGFBP4 null mice have decreased body fat through decreased PPAR-γ expression. In the current study, we assessed whether variation in the IGFBP4 coding region influences body mass index (BMI) in American Indians who are disproportionately affected by obesity. Whole exome sequence data from a population-based sample of 6779 American Indians with longitudinal measures of BMI were used to identify variation in IGFBP4 that associated with BMI. A novel variant that predicts a p.Ser76Thr in IGFBP4 (Thr-allele frequency = 0.02) was identified which associated with the maximum BMI measured during adulthood (BMI 39.8 kg/m2 for Thr-allele homozygotes combined with heterozygotes vs. 36.2 kg/m2 for Ser-allele homozygotes, ß = 6.7% per Thr-allele, p = 8.0 × 10-5, adjusted for age, sex, birth-year and the first five genetic principal components) and the maximum age- and sex-adjusted BMI z-score measured during childhood/adolescence (z-score 0.70 SD for Thr-allele heterozygotes vs. 0.32 SD for Ser-allele homozygotes, ß = 0.37 SD per Thr-allele, p = 8.8 × 10-6). In vitro functional studies showed that IGFBP4 with the Thr-allele (BMI-increasing) had a 55% decrease (p = 0.0007) in FOXO-induced transcriptional activity, reflecting increased activation of the PI3K/AKT pathway mediated through increased IGF signaling. Over-expression and knock-down of IGFBP4 in OP9 cells during differentiation showed that IGFBP4 upregulates adipogenesis through PPARγ, CEBPα, AGPAT2 and SREBP1 expression. We propose that this American Indian specific variant in IGFBP4 affects obesity via an increase of IGF signaling.

Cardiorespiratory Fitness, BMI, Mortality, and Cardiovascular Disease in Adults with Overweight/Obesity and Type 2 Diabetes.

INTRODUCTION: We estimated the effects of cardiorespiratory fitness (CRF) and body mass index (BMI) at baseline on mortality and cardiovascular disease events in people with type 2 diabetes who participated in the Look AHEAD randomized clinical trial. METHODS: Look AHEAD compared effects of an intensive lifestyle intervention with diabetes support and education on cardiovascular disease events in 5145 adults age 45-76 yr with overweight/obesity and type 2 diabetes. In 4773 participants, we performed a secondary analysis of the association of baseline CRF during maximal treadmill test (expressed as metabolic equivalents (METs)) on mortality and cardiovascular disease events during a mean follow-up of 9.2 yr. RESULTS: The mean (SD) CRF was 7.2 (2.0) METs. Adjusted for age, sex, race/ethnicity, BMI, intervention group, and ß-blocker use, all-cause mortality rate was 30% lower per SD greater METs (hazard ratio (HR) = 0.70 (95% confidence interval, 0.60 to 0.81); rate difference (RD), -2.71 deaths/1000 person-years (95% confidence interval, -3.79 to -1.63)). Similarly, an SD greater METs predicted lower cardiovascular disease mortality (HR, 0.45; RD, -1.65 cases/1000 person-years) and a composite cardiovascular outcome (HR, 0.72; RD, -6.38). Effects of METs were homogeneous on the HR scale for most baseline variables and outcomes but heterogeneous for many on the RD scale, with greater RD in subgroups at greater risk of the outcomes. For example, all-cause mortality was lower by 7.6 deaths/1000 person-years per SD greater METs in those with a history of cardiovascular disease at baseline but lower by only 1.6 in those without such history. BMI adjusted for CRF had little or no effect on these outcomes. CONCLUSIONS: Greater CRF is associated with reduced risks of mortality and cardiovascular disease events.

Increased Adiposity and Low Height-for-Age in Early Childhood Are Associated With Later Metabolic Risks in American Indian Children and Adolescents.

BACKGROUND: Growth abnormalities in childhood have been related to later cardiometabolic risks, but little is known about these associations in populations at high risk of type 2 diabetes. OBJECTIVES: We examined the associations of patterns of growth, including weight and height at ages 1-59 months, with cardiometabolic risk factors at ages 5-16 years. METHODS: We linked anthropometric data collected at ages 1-59 months to cardiometabolic data obtained from a longitudinal study in a southwestern American Indian population at high risk of diabetes. Analyses included 701 children with ≥1 follow-up examination at ages 5-16 years. We derived age- and sex-specific weight-for-height z-scores (WHZ) and height-for-age z-scores (HAZ) at ages 1-59 months. We selected the highest observed WHZ and the lowest observed HAZ at ages 1-59 months and analyzed associations of z-scores and categories of WHZ and HAZ with cardiometabolic outcomes at ages 5-16 years. We used linear mixed-effects models to account for repeated measures. RESULTS: Overweight/obesity (WHZ >2) at ages 1-59 months was significantly associated with increased BMI, fasting and 2-hour postload plasma glucose, fasting and 2-hour insulin, triglycerides, systolic blood pressure, diastolic blood pressure, and decreased HDL cholesterol at ages 5-16 years relative to normal weight (WHZ ≤1). For example, at ages 5-9 years, 2-hour glucose was 10.4 mg/dL higher (95% CI: 5.6-15.3 mg/dL) and fasting insulin was 4.29 µU/mL higher (95% CI: 2.96-5.71 µU/mL) in those with overweight/obesity in early childhood. Associations were attenuated and no longer significant when adjusted for concurrent BMI. A low height-for-age (HAZ < -2) at ages 1-59 months was associated with 5.37 mg/dL lower HDL (95% CI: 2.57-8.17 mg/dL) and 27.5 µU/mL higher 2-hour insulin (95% CI: 3.41-57.6 µU/mL) at ages 10-16 years relative to an HAZ ≥0. CONCLUSIONS: In this American Indian population, findings suggest a strong contribution of overweight/obesity in early childhood to cardiometabolic risks in later childhood and adolescence, mediated through persistent overweight/obesity into later ages. Findings also suggest potential adverse effects of low height-for-age, which require confirmation.

Functional variants in cytochrome b5 type A (CYB5A) are enriched in Southwest American Indian individuals and associate with obesity.

OBJECTIVE: This study aimed to identify genetic variants enriched in Southwest American Indian (SWAI) individuals that associate with BMI. METHODS: Whole genome sequencing data (n = 296) were used to identify potentially functional variants that are common in SWAI individuals (minor allele frequency ≥10%) but rare in other ethnic groups (minor allele frequency < 0.1%). Enriched variants were tested for association with BMI in 5,870 SWAI individuals. One variant was studied using a luciferase reporter, and haplotypes that included this variant were analyzed for association with various measures of obesity (n = 917-5,870), 24-hour energy expenditure (24-h EE; n = 419), and skeletal muscle biopsy expression data (n = 207). RESULTS: A 5' untranslated region variant in cytochrome b5 type A (CYB5A), rs548402150, met the enrichment criteria and associated with increased BMI (ß = 2%, p = 0.004). Functionally, rs548402150 decreased luciferase expression by 30% (p = 0.003) and correlated with decreased skeletal muscle CYB5A expression (ß = -0.5 SD, p = 0.0008). Combining rs548402150 with two splicing quantitative trait loci in CYB5A identified a haplotype carried almost exclusively in SWAI individuals that associated with increased BMI (ß = 3%, p = 0.0003) and decreased CYB5A expression, whereas the most common haplotype in all ethnic groups associated with lower BMI and percentage of body fatness, increased 24-h EE, and increased CYB5A expression. CONCLUSIONS: Further studies on the effects of CYB5A on 24-h EE and BMI may provide insights into obesity-related physiology.

A missense variant Arg611Cys in LIPE which encodes hormone sensitive lipase decreases lipolysis and increases risk of type 2 diabetes in American Indians.

AIMS: Hormone sensitive lipase (HSL), encoded by the LIPE gene, is involved in lipolysis. Based on prior animal and human studies, LIPE was analysed as a candidate gene for the development of type 2 diabetes (T2D) in a community-based sample of American Indians. MATERIALS AND METHODS: Whole-exome sequence data from 6782 participants with longitudinal clinical measures were used to identify variation in LIPE. RESULTS: Amongst the 16 missense variants identified, an Arg611Cys variant (rs34052647; Cys-allele frequency = 0.087) significantly associated with T2D (OR [95% CI] = 1.38 [1.17-1.64], p = 0.0002, adjusted for age, sex, birth year, and the first five genetic principal components) and an earlier onset age of T2D (HR = 1.22 [1.09-1.36], p = 0.0005). This variant was further analysed for quantitative traits related to T2D. Amongst non-diabetic American Indians, those with the T2D risk Cys-allele had increased insulin levels during an oral glucose tolerance test (0.07 SD per Cys-allele, p = 0.04) and a mixed meal test (0.08 log10 µU/ml per Cys-allele, p = 0.003), and had increased lipid oxidation rates post-absorptively and during insulin infusion (0.07 mg [kg estimated metabolic body size {EMBS}]-1  min-1 per Cys-allele for both, p = 0.01 and 0.009, respectively), compared to individuals with the non-risk Arg-allele. In vitro functional studies showed that cells expressing the Cys-allele had a 17.2% decrease in lipolysis under isoproterenol stimulation (p = 0.03) and a 21.3% decrease in lipase enzyme activity measured by using p-nitrophenyl butyrate as a substrate (p = 0.04) compared to the Arg-allele. CONCLUSION: The Arg611Cys variant causes a modest impairment in lipolysis, thereby affecting glucose homoeostasis and risk of T2D.